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ABSTRACT 

We make an inventory of the baryonic and gravitating mass in structures ranging from the smallest 
galaxies to rich clusters of galaxies. We find that the fraction of baryons converted to stars reaches 
a maximum between M500 = 10 12 and 10 13 M Q , suggesting that star formation is most efficient in 
bright galaxies in groups. The fraction of baryons detected in all forms deviates monotonically from 
the cosmic baryon fraction as a function of mass. On the largest scales of clusters, most of the expected 
baryons are detected, while in the smallest dwarf galaxies, fewer than 1% are detected. Where these 
missing baryons reside is unclear. 

Subject headings: cosmological parameters — dark matter — galaxies: clusters: general — galaxies: 
dwarf — galaxies: irregular — galaxies: spiral 
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1. INTRODUCTION 

The early universe was a highly uniform and homo- 
geneous mix of dark and baryonic matter with baryon 
fraction f b = 0.17 ±0.01 ([Komatsu et al.l [2009h . If 
this primordial mix persists as individual gravitation- 
ally bound structures emerge during the course of cos- 
mic evolution, then the baryonic mass of any given ob- 
ject would be Mb = fbM to t- Indeed, combining this with 
the ba ryon density constra int from big bang nucleosyn- 
thesis (( Walker et al.l 119911 ) provides one important ar- 
gument that the density parameter is less tha n unity 
([White et afl 119931 : lOstriker &; Steinhardtl 119951) . Here 
we reverse the logic and ask what fraction of the expected 
baryons are actually detected. 

We present an inventory of the baryonic and gravitat- 
ing masses of cosmic structures spanning a dozen decades 
in detected baryonic mass. The mass of baryons known 
in each system correlates well with the total mass, but 
not as a simple proportion. The implication is that most 
of the baryons associated with individual dark matter ha- 
los are n ow missing. While some are in the intergalactic 
medium ([Danforth &; Shul3 l2005). a complete accounting 
of where the baryons now reside, and how they relate to 
their parent structures, remains wanting. 

2. MASS INVENTORY 

We divide our inventory into two broad categories of 
gravitationally bound systems: those supported by ro- 
tation, and those supported by random motions. This 
distinction is important to the way in which we infer the 
gravitating mass, as discussed for each type of system be- 
low. The detected baryonic mass is the sum of observed 
stellar and gaseous components: M& = M* +M g . In gen- 
eral, the gas mass is more precisely known as the physics 
of the emission mechanism is better understood, while 
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the stellar mass requires an estimate of the mass-to-light 
ratio of a composite stellar population in order to convert 
observed luminosity to mass. Over the many decades in 
mass considered here, a factor of ~ 2 uncertainty in the 
the stellar mass-to-light ratio is a minor concern. 

2.1. Rotationally Supported Systems 

One of the most important lines of evidence for 
dark matter is the observation that the rotation curves 
of spiral galaxies tend to bec ome asymptotically flat 
(|Rubin et alj|1978t lBosmal[l98l[ ) at large radii where the 
baryonic mass alone should result in a Keplerian decline. 
This outermost circular velocity, V c , we take to be rep- 
resentative of the gravitating mass. Among rotating sys- 
tems, we distinguish between those dominated by stellar 
mass (typically early type spirals) and those dominated 
by gas mass (typically irregular and late type spirals). 

2.1.1. Star Dominated Spiral Galaxies 

Bright spiral galaxies have the majority of their de- 
tected baryonic mass in the form of stars. We utilize 
the ro tation velocities and mass estimates of IMcGaughl 
( 2005), selecting only those galaxies with M* > M g for 
Q = 1. The r esults are in g ood agreement with inde- 
pendent work (de Blo k et all [20081 ). The sum of bary- 
onic mass is completed by adding the cold gas mass, 
corrected by a factor of 1.4 to account for the presence 
of helium and heavier elements. Dust and hot gas do 
not contribute significantly to the total (|Bregmanll2007t 
lAnderson fc B regman 200^). 

2.1.2. Gas Dominated Galaxies 

Many dim, late type disk galaxies have more of their 
baryonic mass in the form of atomic gas than in stars: 
M* < M g . These arc particularly interesting for our 
purposes here because their baryonic mass is insensi- 
tive to the choice of stellar mass estimator. For these 
galaxies, Mb follows almost directly from the observed 
21 cm flux and the physics of the hydroge n hyperfine 
transition. Two recent indepen dent studies (|Stark et al.l 
l2009t iTrachternach et aT1l2009f) provide for the first time 
a large combined sample of such galaxies. 

Rotationally supported disk galaxies are known to 
obey a relation between baryonic mass and rotation 



2 



M 500 ( M ©) 

7 8 9 10 11 12 13 14 15 

2 10 10 10 10 10 10 10 10 10 
o r — 1 r-n-*i "i — i i i h — i i i 1 1 1 % Yi 




03 

O L i i i i 



10 1 10 Z 10 3 

V c (km s" 1 ) 

Figure 1. The relation between baryonic mass and rotation ve- 
locity. The sum of detected baryonic mass is plotted against the 
circular velocity of gravitationally bound extragalactic systems. 
Round symbols represent rotationally supported disks while square 
symbols represent pressure supported systems. Larger symbols 
correspond to systems whose baryonic mass is dominated by gas 
and smaller symbols those domina ted by stars. D ark blue cir- 
cles are for star d ominated spirals (McGaugh 2005). Light blue 
( Stark ct al. 2009]) and green ( Trachtcrnach ct al. 2009) circles rep- 
resent recent work on gas dominat ed disks. Red squa res repre- 
sent Local Group dwarf satellites ([Walker ct al. 20 0lj). Purple 
squar es represent the mean of many galaxy clusters HGiodini et al.1 
2009). These pressure supported systems fall close to, but system- 
atically below the Baryonic Tully-Fisher relation defined by the 
disks (dashed line). The mass enclosed by an over-density 500 
times the critical density is shown by the upper abscissa assum- 
ing M500 = (1.5 x 10 5 km -3 s 3 Mg )V C 3 (see discussion in text). 
If structures possessed baryons in the cosmic fraction (/), = 0.17), 
they would fall along the solid line. 

velocity known as the Baryonic Tully-Fisher relation 
(|McGaughll200l . This is shown in Fig. [TJ giving empha- 
sis to the gas dominated galaxies (larger circles) which 
extend the relation beyond that previously known for 
bright spirals. The y also provide a physics-based calibra- 
tion o f the relation ([Stark et al.ll2009llTrachternach et alJ 
2009) . This is shown as the dashed line i n Fig, ffl The in - 
trinsic relation appears to be very tight (|Verheiienl r2001'). 
with observational uncertainty (due mostly to the un- 
certainty in the distances to the individual galaxies) ac- 
counting for essentially all of the scatter. 

The Baryonic Tully-Fisher relation for rotating disks 
shown in Fig. [T] extends over five decades in baryonic 
mass, a considerable improvement over the two decades 
typically considered. Before considering how the bary- 
onic mass relates to the total mass, we first investi- 



Table 1 

Broken Power Law Fit 



Scale 


Range in V 


X 


A 


Dwarf 


< 20 


8.2 


-3.88 


Disk 


20 — 350 


4.0 


1.65 


Cluster 


> 350 


3.2 


3.69 



Note. — log Mb = x log V c + A 



gate whether this relation can be extended still fur- 
ther. Rotationally supported systems cover the range 
20 < V c < 300 kms" 1 , but there are both smaller and 
larger pressure supported systems. 

2.2. Pressure Supported Systems 
2.2.1. Elliptical Galaxies 

Observations of giant elliptical galaxies typically do 
not exten d far enough radially to identify the equiva- 
lent of V c (Cappcll ari et al.H2007f). Gravita tional lensing 
(jHoekstra et al.l 120051 : IGavazzi et all 120071 ) provides im- 
portant constraints on the total mass, but the conversion 
to the equivalent circular velocity is sensitive to the as- 
sumed model. It is also subject to degeneracy between 
stellar and dark mass. Given these uncertainties, and 
that these data largely overlap with those for spirals, we 
do not consider giant ellipticals furthe r. We do note that 
the es timates of Hockstr a et al.1 ([2005) and Gavazz i~et al.1 
(2007) are consistent with our results to the same degree 
that they are consistent with each other. 

2.2.2. Local Group Dwarfs 

Recent years have witnessed an explosion in the dis- 
covery of small satellite galaxies of both the Milky Way 
(Bclo kurov et all 120071 ) and M31 (|McConnachie et all 
l2008f ). The proximity of these quasi-spherical satellite 
galaxies makes it possible to identify isolated stellar sys- 
tems much smaller than known elsewhere, and to mea- 
sure their inte rnal kinematics wi t h velocities from in- 
divid ual stars (| Walker et al.l [20071 [20091 : iSimon fc Gehal 
|2007[ ). We assume isotropic internal orbits to relate the 
observed line-of sight velocity dispersion a to the circu- 
lar velocity: V c = V3<7 ([Wolf et al.ll2009h . The baryonic 
masses of these systems are dominated by stars. To con- 
vert light to mass, we adopt the mean s tellar mass-to- 
light ratio estimated (|Mateo et al.l I199H ) from resolved 
stellar population studies: M^/Ly — 1.3 M /Lq. 

2.2.3. Clusters of Galaxies 

Recent work (|Giodini et al.|[2009[ ) on clusters provides 
a large, homogeneous set of data with stellar, gas, and 
gravitating masses averaged over many individual clus- 
ters. Most of the baryonic mass in clusters is in hot, X- 
ray emitting gas. This provides a good measure of both 
the gas mass and the gravitati ng mass. Stellar ma sses 
are based on a scaling relation (|Giodini et alJ [20091 ) de- 
rived from if-band data. The result here is not sensitive 
to the precise relation as the X-ray gas dominates the 
baryon budget. 

It is conventional to refer to structures by the den- 
sity contrast they represent with respect to the criti- 
cal density of the universe. The mass enclosed within 
a radius encompassing the over-density A is A^a = 
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Table 2 

Binned Data 



System <V C > ay <M b > a b <M*> cr, <Msoo> fd <* f f* CT * Ref. 
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References. 1. Clu sters: IGiodini et "all J2009l>; 2. Spirals: IMcGaughl (t2005tl : Gas dominated galaxies: 3. IStark et all ( f2009^ and 4. 
ITrachternach et al.l [120091) : 5. Local Group dwarfs: IWalker et alj 12009) . 

Note. — Velocities and masses are logarithmic with units of km s -1 and solar masses, respectively. 



(47r/3) Ap cr i t R&- With the definition of critical den- 
sity and circular v elocity, M A = (A/2 )- 1 / 2 (GH )~ 1 V^. 
The cluster data (|Giodini et al.l [2009h are referenced to 
A = 500, so for H a = 72 kms" 1 Mpc -1 (jFreedman et al.l 
[200l . M 500 = (20 4 5 5 2 km -3 s 3 Mq 1 )^. To relate the 
circular velocity V c observed in spirals to V500 at larger 
radii, we note that rotation curves are approximately 
flat, so we expect V c = fvVsoo with fy « 1. If wc 
associate the observed V r with t he peak velocity of an 
NFW halos (jNavarro et al.lfl99l . then 1.0 < fv < 1-3 
over the relevant range of masses. This appears consis- 
tent with the Milky W ay (|Flvnn et all 120061 : Kue et all 
[2001 IMcGaughl [2001 if f v = 1.1, which we adopt for 
specificity. This is also cons istent with the independent 
estimate of ILau et al.l ([2009) . Note that a 20% change in 
fv corresponds to shifting the cluster data by the width 
of a data point in Fig. [TJ 

The choice of the radius -R500 is driven by the data, 
and is inevitably somewhat arbitrary. A preferable choice 
might be th e virial radius, which oc curs around A « 100 
in ACDM (jBrvan fc Normanl fl99ct ) . The data do not 
constrain this, nor is there anything genuinely special 
about the virial radius as the halo profile is expected 
to merge smoothly with the cosmic background. So we 
reference our masses to i?sooi and caution that other 
plausible choices would give apparently conflicting re- 
sults merely because of the differing convention. We also 
note that cluster baryon fractions tend to be rising at the 
last measured point, which might help to reconcile them 
with the cosmic baryon fraction. In contrast, choosing a 
larger (e.g., virial) reference radius for individual galax- 
ies would only include more dark mass without any ad- 
ditional baryons, driving their inferred baryon fractions 
to smaller values. 



2.3. Generalized Baryonic Tully-Fisher Relation 

Fig.[T]represents a generalization of the Baryonic Tully- 
Fisher relation to pressure supported systems, and exten- 
sion to higher and lower mass. The data can be described 
by a broken power law, with a different slope and nor- 
malization at the dwarf, disk galaxy, and cluster scale 
(TableQ]). The variation between these regimes is a small 
effect compared to the deviation from a constant baryon 
fraction (solid line in Fig. [T]) . The gravitating mass scales 
as M ~ V 3 , but the observed relation is steeper at all 
scales. Consequently, the detected baryon mass deviates 
systematically from the expectation of a constant cosmic 
baryon fraction. 

2.4. Binned Data 

The data of IGiodini et alJ (|2009l ) represent the bin- 
ning of many individual clusters, while the remaining 
data in Fig. [T] represent individual galaxies. To place 
all data on the same footing, we also bin the galaxy 
data. We bin in intervals in baryonic mass containing 
roughly equal numbers of galaxies. We maintain the dis- 
tinction between galaxies of different types, and for the 
gas dominat ed galaxies, we ma intain separation between 
the da ta of IStark et "ail (|2009f ) and ITrachternach et all 
(2009) . The binned data are reported in Table [2] and 
plotted in Fig. [2j Table [2] gives the mean and stan- 
dard deviation in the mean of each logarithmic bin in 
circular velocity, baryonic mass, stellar mass, and grav- 
itating mass within A = 500. It also gives the detected 
baryon fraction fd = Mb/ifbMsoo), which is the ratio of 
known baryons to the amount expected from the cosmic 
baryon fraction. The fraction of these nominally avail- 
able baryons that have been converted into stars is given 
by the stellar fraction /* = M^Kf^M^o). 
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Figure 2. The fraction of the expected baryons that are detected [left: fd = Mb/(ff,Msoo)] and the fraction converted to stars [right: 
/* = M*/(/(,Msoo)]' Symbols have the same meaning as in Fig.[T]but with the binning of Table[2] The detected baryon fraction increases 
monotonically with mass while the stellar fraction peaks between M500 = 10 12 and 10 13 Mq. 



3. DISCUSSION 

The stellar fraction reaches a maximum between 
Af 500 = 10 12 and 10 13 M Q (Fig.©. This is broadly con- 
sistent with previo us results based on counting statistics 
(|Yang et alj I2008T ) . However, there may be some offset 
in the peak mass scale relating to the long standing di- 
chotomy between Tully-Fisher and luminosity function 
based normalizations of the halo mass function. Notably 
the efficiency of star formation appears to i ncrease mono- 
tonic ally with mass for individual galaxies (jBaldrv et alJ 
2008). After the transition to cluster halos containing 
many galaxies, the efficiency declines again. 

Perhaps the most striking aspect of Fig. [2] is that the 
fraction of detected baryons falls short of the cosmic frac- 
tion at all scales. In no system is it unity as would 
be expected if we had a complete accounting of all the 
baryons associated with a given bound system. Where 
are all these missing baryons? 

An obvious possibility is that the missing baryons are 
present, and simply inhabit their dark matter halos in 
some undetected form. Indeed, one might expect that 
not all baryons would have time to cool into the observed 
cold gas and stellar component of galaxies. In this case, 
many baryons might remain mixed in with the dark halo. 

The notion that we are simply not seeing many or 
even most of the baryons in individual galaxies is pro- 
foundly unsatisfactory. Direct searches for hot halo gas 
have turned up nothing substantial: halo baryon reser- 
voirs fail to explain the observed deficit by two orders o f 
magnitude ([Bregman 2007t lAnderson fc B regman 2009) . 
In clusters, the hot gas is detected, and constitutes the 
majority of the baryons. Clusters fall sh ort of the cosmic 
baryon fraction by a m odest amount ((McCarthy et al.l 
l2007|: | Giodini et al]|2009f) which might be readily explica- 
ble (jCrain et al.ll2007T ). However, the fraction of missing 
baryons is highly significant on the scales of individual 
galaxies. In dwarfs, fewer than 1% of the baryons ex- 
pected from the cosmic fraction are detected. 

The detected baryon fraction varies systematically 
with scale. It is a matter of taste whether one chooses 



to describe this scale as one of circular velocity, mass, 
or potential well depth. It is tempting to attribute 
this correlation to feedback processes bein g more effec- 
tive i n objects with smaller potential wells (jDekel fc Silkl 
1986). H owever the details are heinously complicated 
(jMaver fc Moorei 120041 : IKeres et alj I2009D and not well 
understood. We would naively expect feedback to be 
a messy process resulting in lots of scatter in any cor- 
relations that might result. Instead, the observed rela- 
tion is remarkably tight. In principle, the entire range 
< fd < 1 is accessible at each mass, yet only a very 
particular value is observed. Moreover, the current po- 
tential well is the result of the hierarchical assembly of 
many smaller building blocks. Left to itself, a small dark 
matter halo will have a small fd- If incorporated into a 
larger halo, fd goes up. How does each building block 
know to bring the right amount of baryons to the final 
halo? 

We do not see a satisfactory solution to this missing 
baryon problem at present. Considerable work remains 
to be done to obtain a complete understanding of the 
universe and its contents. 
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